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Abstract. A 10-night asteroseismic observation programme has been conducted in January 2004 with the spec- 
trometer HARPS at the ESO 3.6-m telescope. The selected target, the 6th magnitude F5V star HD 49933, was 
chosen among the prime candidates of COROT, the European space mission dedicated to characterize stellar os- 
cillations mode with high precision photometry measurements. This star shows important line profiles variations, 
indicating a surprisingly high activity with respect to its low rotation rate. However, with the help of tools devel- 
oped for disentangling the signatures of activity and oscillations, we are able to observe its oscillation spectrum in 
the frequency range [1.2, 2.2 mHz]. We measure the large separation (88.7±0.4 //Hz) and the maximum amplitude 
(around 0.4 ± 0.1 ms~^ rms), respectively in agreement and marginal agreement with the predicted values. 
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' 1. Introduction 

F type stars lack seismic observations owing to the diffi- 
culty, due to the poor number of their spectral lines, to 
, get a high accuracy seismic signal with ground-based spec- 
■ trometric measurements. In this study, we present an as- 
, teroseismic observational programme conducted with the 
' new generation spectrometer HARPS, installed at the ESO 
3.6-m telescope in La Silla (Pepe et al. I2002|) on a F type 
' star which is a candidate target (HD 49933) of the space 
COROT mission (Baghn et al. I2002|) . The unique capabili- 
ties of the spectrograph HARPS made possible to measure, 
to identify and to characterize the solar-like oscillations of 
such a 6th magnitude F star. 

In this letter, we present first the observations 
(Sect. 2), then describe the long term radial velocity pat- 
tern we observed (Sect. 3). The method we used in or- 
der to correctly extract the radial velocities is exposed in 
Sect. 4. The seismic signature is confirmed and discussed 
in Sect. 5. Section 6 is devoted to conclusions. 
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* Based on observations obtained with the HARPS echelle 
spectrometer mounted on the 3.6-m telescope at ESO-La Silla 
Observatory (Chile) 



2. Observations 

HD 49933 (HR 2530, ADS 5505A) was selected among 
the COROT prime targets. This F5V star lies in a re- 
gion of the HR diagram not yet explored by previous seis- 
mic observations. HD 49933 has an effective temperature 
of 6500 ± 100 K with an estimated absolute magnitude 
Mv = 3.39 ± 0.10. This star, with an iron abundance of 
—0.37 dex (Solano et al. l2UU4|l . is sligthly metal poor com- 
pared to the Sun and to Procyon. From observations of 
solar-like oscillations in Doppler velocity, Samadi et al. 
(2004| find that the maximum in the oscillation velocity, 
Vmax, scales approximatively as (L/M)°-^^ where L and M 
are the luminosity and the mass of the star respectively. 
Such a scaling law then predicts Knax — 60 ± 0.5 cms~^ 
rms for HD 49933. Spectrometric seismic observations be- 
ing limited by the stellar line width, we had to select a 
slow rotator: HD 49933 had an estimated v sini about 10.9 
kms-i (Solano et al. 2004), changed to 9.5±0.3 kms^^ 
with our more precise measurements. 

About 57 hours observations were recorded between 
January 15 and 25, representing 1351 measurements. The 
exposure time was 120 s, giving one measurement each 
153 s when including overhead, and a Nyquist frequency 
about 3.3 mHz above the predicted cutoff frequency at 
2.2 mHz. The mean SNR in the spectrum (at 550 nm and 
with a spectral bin of 0.85 kms~^) typically varies from 60 
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Fig. 1. Time series of the radial velocity of HD 49933 
measured with the pipehne reduction of HARPS. The aver- 
age line bisector amplitudes, muhiplied by the sign of the 
mean slope of the bisectors, are superimposed (in grey) 
and translated for display purposes. The full signal (up- 
per window) shows large variations, with a period about 
7.9 days. The high frequency component (lower window) 
was corrected from the low frequency fit (thin solid line in 
the upper window) and also from the residual offsets and 
slopes. 



at high airmass to about 120. This SNR corresponds ex- 
actly to the contribution of the photon noise, which varies 
rapidly with the airmass, as the input flux in the entrance 
fiber. The mean observation duration did not exceed 7.5 
hours per night at this summer period. 

3. Long term radial velocity pattern, line 
asymmetries and activity 

The radial velocity curve of HD 49933, determined with 
the online pipeline reduction of harps, shows large am- 
plitude variations with a characteristic time scale of a few 
days (Fig. Further measurements with the CORALIE 
spectrometer on the Swiss 1.2-m telescope (Queloz et al. 
I2000|l revealed that this long term radial velocity pattern 
is not stable over periods longer than 10 days, and cannot 
be attributed to one or several orbiting giant planets. 

Suspecting that long term variations may be the sig- 
nature of stellar activity, we examined the spectral line bi- 
sectors in our spectra. We first constructed average photo- 
spheric profiles, combining the profiles of the many lines in 
the HARPS spectra using the Least-Square Deconvolution 
(LSD) method described in Donati et al. H1997fl . the line 
pattern being constructed using a Kurucz model with 
Teff = 6500 K, log 3=4.5, and [Fe/H]=-0.3. These param- 
eters represent the closest choice to those of HD 49933 in 
the model grids we have at our disposal. We have checked 
that using adjacent models in the grid does not modify 
the results. Typical signal-to-noise ratios of 600 to 700 are 
obtained per velocity bin of 0.2 kms^^ in the LSD average 
profiles. 
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Fig. 2. Average line bisector amplitude of HD 49933. The 
large variations identified in Fig.^are related to the stellar 
activity modulated by the rotation of the star. 



The bisectors of the resulting average line profiles were 
then constructed as the midpoints of horizontal segments 
bounded by the sides of the profiles. We find that the line 
bisectors of HD 49333 most often exhibit the classical C- 
shape of cool stars, although with a larger distortion in the 
lower part than in the upper part of the line, in contrast 
to most other observed cool star bisectors (Gray 1992). 
They have amplitudes up to about 1 kms^"'^, measured as 
the velocity span covered by the bisector, and they are 
highly variable. The bisector variation is well visible on 
Fig. 121 with the 8-day periodicity suggested by the radial 
velocity curve in Fig. ^ We conclude that this long-term 
radial velocity variation can be entirely attributed to line 
profile distortions, which are most likely related to stellar 
surface activity. The 8-day periodicity may be identified 
as rotational modulation of the surface pattern giving rise 
to the variations. One probable interpretation is that the 
line distortions are due to the presence of one or several 
cool spots at the stellar surface, whose signature crosses 
the line profiles from blue to red during the star rotation, 
thus affecting the line bisectors. This interpretation is fully 
compatible with the apparent double-periodicity of the 
signal. The bumps produced by stellar spots are stronger 
near the bottom of the lines than in the line wings, which 
may explain the stronger variations of the line bisectors 
near their lower end. 

Line bisector variations have already been reported 
in several cases, mostly for G and K type stars (Lopez- 
Santiago et al. EU03. and references therein), the bisector 
modulation being interpreted as the signature of one or 
several cool starspots at the stellar surface. The observed 
bisector modulation amplitudes are usually small. In com- 
parison to these results, the larger amplitude (1 kms^^) 
for the bisector modulation of HD 49933 is somewhat 
surprising for such a slow rotator (wsini = 9.5 kms^^. 
Plot — 8 days). The detailed interpretation of these re- 
sults and the modelling of the stellar surface structures 
will be proposed in a further paper. 
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4. Disentangling Doppler shifts and line profile 
variations 

The signature of activity occuring at low frequency should 
not affect the seismic signal at higher frequency. However, 
the large distortions of the lines related to activity may 
mimic rapidly varying tiny Doppler shifts. Therefore, any 
method for extracting the radial velocity must be able to 
distinguish between a line distortion and a pure line shift. 
The numerical cross-correlation technique (Baranne et al. 

of the reduction pipeline of harps uses as an input 
the theoretical position of the stellar lines. The mean line 
profile is then reconstructed by interpolation of the ob- 
served lines corresponding to the theoretical positions. A 
Gaussian fit gives the Doppler shift, which is affected by 
the line profile. Therefore, deformations of the line pro- 
file may give a spurious Doppler-like signal. The optimum 
method fConnes ll985[ Bouchy et al. I2001|l is in principle 
independent of the line profile. Any spectral line is com- 
pared to itself at different times of observations and the 
wavelength shift is translated into a Doppler shift accord- 
ing to a first order expansion. The efficiency of the method 
requires a reference spectrum with a high SNR. With the 
present set of observations, the method proves to be effi- 
cient in practive with a reference spectrum having a SNR 
greater than 500, resulting of the addition of more than 
30 spectra. Therefore, the line profiles of the reference and 
of the spectra do not coincide, so that Doppler artefacts 
may be induced by lines distortions. Finally, both reduc- 
tion methods give time series with similar velocity noise 
of about 3.5 ms~^ per integration time (Fig.^. However, 
we cannot be sure that the signal in the Fourier spectrum 
plotted in Fig. |31is not appreciably contaminated by the 
stellar activity. 

In order to disentangle the distortion and Doppler shift 
as precisely as possible, we have analysed the LSD profile 
by decomposing it in two terms. The mean line profile 
I{t) being considered as a function of the velocity v, the 
intensities I(t,v) and I{t,—v) explore it symmetrically, 
and give the integrated quantities: 

V 

I+,^{t) = Y.w{v) [l{t,v)±I{t,-v)] 



The normalized weighting function w{v) is calculated from 
a Gaussian fit of /, with w proportional to the derivative 
of this fit, this enhancing the role of the line wings com- 
paratively to the continuum and the line center. The ve- 
locity V is high enough for covering the whole profile. The 
spectrum of the symmetric component /_|_ (t) shows no ex- 
cess power in the range [1.2, 2.2 mHz] as seen on Fig. 0] 
contrary to the antisymmetric signal /_ (t) , which is very 
similar to the spectrum of Fig. O /+ is sensitive to the line 
broadening and distortion, while /_ is not, but sensitive to 
a line shift. This treatment is then able to extract the in- 
formation of a pure Doppler signal, and confirms that the 
observed Fourier spectrum is not affected by distortions. 
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Fig. 3. Fourier spectrum of the velocity time series cal- 
culated with the optimum method. The window function 
is shown in the inset. No high pass filtering was applied, 
except for the removal of the mean velocity slope of each 
night. 
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Fig. 4. Fourier spectrum of the time series /+ (upper 
curve) and /_ (lower curve), with a signal translated in 
velocity. The spectrum of /+ , representative of the defor- 
mation of the mean line, does not show any excess power 
in the frequency range [1, 2.2 mHz], contrary to the one 
of very similar to Fig.|31 



5. The seismic signature 

The Fourier spectrum shows an excess power in the fre- 
quency range [1.2, 2.2 mHz], peaking at 1.7 mHz near 
the predicted value of 1.3 mHz derived from the scal- 
ing law Bedding & Kjeldsen H2003|l . The comb analysis 
(Mosser et al. [iy981 exhibits a clear signature around 
Aiy = 88.7 ± 0.4 /iHz (Fig. 0, representative of the seis- 
mic large separation. We computed evolutionary models 
for HD 49933, using the CESAM code fMorei nWTj l. We 
considered the global parameters of HD 49933 (6500 K, 
[Fe/H] = —0.37, My — 3.39) and classical physical op- 
tions for this kind of objects. Heavy element abundance 
Z= 7 10"'^ is deduced from [Fe/H] and we adopted for he- 
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Fig. 5. Comb response of the spectrum of /_ between 1.2 
and 2.2 mHz. The signal around 88.7 //Hz is the signature 
of low degree solar-like pressure modes. 

hum Y=0.255 with Yprimordial=0.240 at Zprimordial = 0,with 

Z/X© = 0.0245 and an enrichment galactic law dY/dZ— 2. 
The low metallicity induces a small mass for this range of 
temperature, so that the preferred mass is around 1.15 
Mq, with an age between 3.3 and 4.0 Gyr. For all com- 
puted models, the associated seismic large separation is 
found to vary between 75 and 88 /iHz. The 88.7 /iHz spac- 
ing found in the data is thus consistent with this range of 
expected values, and is in favor of a young age. 

In order to estimate the maximum oscillation ampli- 
tude, we have constructed synthetic time series, based 
on a theoretical low degree p-modes eigenfrequency pat- 
tern. The modes lifetimes are estimated from the eigenfre- 
quency widths, between 1 and 4 /iHz FWHM (Houdek et 
al. 11999(1 . The maximum amplitudes are assumed to follow 
a Gaussian distribution in frequency. The synthetic power 
spectra are then calculated using the model of a stochasti- 
cally excited, damped harmonic oscillator (Anderson et al. 
11990(1 . We assume that the spectrum is composed of this 
synthetic signal and a white noise. The maximum ampli- 
tude was then determined in order to obtain comparable 
energy per frequency bin in the synthetic and observed 
spectra. A Monte-Carlo approach finally shows that the 
best agreement is for a signal maximum amplitude about 
0.4±0.1 ms~^ rms, with a Gaussian envelope centered at 
1.7 mHz and with a 1.0 niHz FWHM. This estimate is not 
much sensitive to the assumed eigenfrequency widths: an 
uncertainty of ±50% of their values adds an uncertainty 
smaller than 15%. The agreement with the predicted scal- 
ing (60±5 cms^^ rms) appears to be marginal. 

6. Conclusion 

We have been able to measure with HARPS the p-modes os- 
cillation amplitude and the large separation of HD 49933. 
This confirms that ground-based asteroseismologic obser- 
vations with HARPS provide excellent solar-like oscillations 
studies, even for a faint (my — 5.7) F type star with 



a tisini near 10 kms^^ and exhibiting activity. For a 
star such as HD 49933, ground-based spectrometric mea- 
surements are less efficient than for cooler types stars, 
with more and narrower lines. However, we have demon- 
strated that spectrometric measurements are the appro- 
priate tools for disentangling the signatures of oscillations 
and activity. This problem is particularly important for 
intermediate mass stars, as most of the COROT targets, 
that contain information on the extension of the core and 
on the angular momentum transfer in the early phases of 
their lives. The question of contamination of the seismic 
signal by activity in the p-mode frequency range has been 
raised by Matthews 1(2004(1 as a possible cause of the non 
detection of oscillations on Procyon with MOST. Our ob- 
servations allow an unambiguous distinction between the 
stellar activity and the seismic signal. 
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